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Abstract 



o 
o 

Q\ ■ A dynamic equation for velocity structure functions S n (r) =< (u(x) — u(x')) n > 



in strong turbulence is derived in the one-loop (eddy viscosity) approximation. This 
homogemeous differential equation yields scaling exponents £ n in the relations S n (r) cx 



', which are in a very good agreement with experimental data. 

O 

Kolmogorov's (1941) relation for the third-order structure function which is in the limit of 

•i-H . 

zero viscosity u — > 0, 



S 3 =< (u(2)-u(l)) 3 >=-Ur 

5 



where 8 = u(diUj) 2 is the mean dissipation rate of turbulent kinetic energy, u(i) = m(xi) is 
the value of the rc-component of the velocity field at the point x, and r = n • (x2 — Xi) is 
the value of the displacement along the x-axis (n is the unit vector in the x-direction) [1]. 
Applying dimensional considerations to this dynamic relation, Kolmogorov also made the 
prediction: 



. 2 2 



S 2 =< ( u (2) - u{l)Y >oc £3 r s (2) 
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leading to the celebrated Kolmogorov energy spectrum: 



E(k) = C K £h~t 



(3) 



with £ = 5/3. Kolmogorov went even further by generalizing (2) to the structure function 
S n (r) of arbitrary order n: 



Experimental investigations have supported the relation (2) with good accuracy, giving 
for the exponent of the second-order structure functions £ 2 « 1.66 — 1.75. At the same 
time substantial deviations of £ n from the K41 values £ n = n/3 from n > 3 have also been 
observed. It is interesting that no data on the scaling exponents of the structure functions 
S n with n < 1 have yet been reported. 

Ever since Kolmogorov's work, dimensional considerations were the most dominant way 
of evaluating the scaling exponents. This is why until recently no progress in understanding 
the intermittency of strong turbulence was achieved. In 1994 R. H. Kraichnan, considering 
the problem of a passive scalar advection in a rapidly- changing- in- time random velocity 
field, realized that anomalies, coming from the dissipation terms in the equation of motion, 
result in a homogeneous differential equation for S n (r), leading to an algebraic form of the 
scalar structure functions S n (r) oc r^ n with scaling exponents £ n determined by the numerical 
values of the coefficients in the equation [2]. This is why the scaling exponents £ n cannot 
be obtained on the basis of dimensional considerations. Later, groundbreaking works by 
Gawedzkii and Kupiainen [3 ], Chertkov et. al. [4 ] proved that in the vicinity of the 
gaussian limits the scaling of scalar structure functions is determined by the zero modes 
of the homogeneous differential equations for S n (r), which for this problem can be written 
explicitly [5 ], [ 6 ]. The results of Refs. [ 3]- [4 ] have been confirmed in ref. [7 ] using a 
different approach. At about the same time Shraiman and Siggia [8 ], considering a model 
of a passive scalar based on the concept of eddy diffusivity, showed that zero modes play the 
most important part in determination of the scaling exponents of structure functions. 



S n (r) =< (u(2) - «(!))" >=A n {Sr)^ 



(4) 
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Since Richardson's (1926) work, effective transport coefficients (eddy viscosity and eddy 
diffusivity) , though not rigorously justified, have been used for the quantitative description 
of transport phenomena in strongly turbulent engineering flows. It is safe to state that by 
now these concepts have evolved into an engineering tool widely used for design purposes 
throughout mechanical engineering. The idea behind the method is simple [9], [10]. Consider 
the Navier-Stokes equations for an incompressible fluid: 

u it + u • Vui = Si- ViP + v \7 2 Ui (5) 
V ■ u = 

where S(x, t) (V • S = 0) is a large-scale source function. Using the incompressibility con- 
dition the pressure p can be expressed in terms of the velocity field. It is useful to write the 
equation for the Fourier-transform u(k, uj): 

-iuui(k, u) + ^Pi mn (k) J dqdQu m (<i, u)u n (k - q, u - Q) = Si(k, u) - v k 2 ui(k, u) 
where the projection operator P im „(k) is: 

-Pimn(k) = A; m -P; ra (k) + k n Pi m (\t) 

and Pij = 5ij — ^f. In the eddy viscosity approximation the effects of the small scales 
(q >> k) on the modes u(k) is represented in terms of the eddy viscosity [10 ] (we neglect 
the small eddy noise contribution): 

-iujui(k, u) + ^Pi mn (k) J dqdilu m (q, u)u n (k - — w Si -Tk 2+a ui(k, uj) - u k 2 ui(k, u) 

(6) 

where 



2 + 6 m 
a = 7T- ( 7 ) 
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The integration is carried out over the interval: < q < k; — oo < uj < oo. This 
means that this equation describes the velocity field averaged over small-scale fluctuations 
with q > k. This result must be considered as a general model for the velocity field ui(k,u) 
where the effects of the small- scale fluctuations are accounted for by the eddy viscosity 
with v(k) oc Yk a . In principle, this equation must be solved subject to initial and boundary 
conditions for each mode w(k, t). A similar model has recently been used for evaluation of 
the anomalous scaling in the problem of a passive scalar advected by a random velocity field 
by Shraiman and Siggia [8 ]. The role of the large scale motions in the dynamics of the 
small-scale velocity fluctuations, described by the non-linear contribution to (6), is two-fold. 
First, the large-scale structures kinematically transfer ( sweep ) the small-scale fluctuations. 
This process does not lead to the energy redistribution. Secondly, the large- scale motions 
serve as an energy source for the small-scale dynamics which can be accounted for by the 
random forcing function f which is a complex functional of the velocity field. Eventually, we 
will be interested in the equation of motion for the structure functions of velocity differences 
Au = u(2) — u(l), for which sweeping is not too important. This allows us to write (6) as a 
simple Langevin-like equation for u(k): 

-ium(k, u)KSi + fi- Fk 2+a Ul {k, uj) - u k 2 Ul (k, uj) (8) 

In this approximation the equation for the structure function S n (r) can be written readily 
in the three-dimensional case: 

^-^ = nD 

where 

D =< (AS + Af + i/ (Vlu(2) - V?«(l)))(A«) n " 1 > 

We are interested in the behaviour of the structure functions for small values r/L « 1. 
Since the source S is assumed to act at large scales only we neglect it in what follows. As 
in the problem of a passive scalar, the derivation of the D-term is a difficult task. Here 
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(9) 



the problem is even harder since we do not know much about the effective energy source f . 
It is clear, however, that the eddy viscosity approximation can be accurate only for S n (r) 
with the relatively small n > 0. Thus, the values of S n are dominated by the part of the 
probability density P(Au,r) where Aw (Au) rms . We assume: 

D = f(r)S n (r) (10) 

with /(r) independent on n. A more detailed argument leading to (10) will be presented 
below where it will be shown that (10) is consistent with the eddy viscosity approximation. 
The expression (10) is not dissimilar to the anzatz introduced by Kraichnan in his theory 
of a passive scalar [2]. The function f(r) is fixed by the relation (1). Introducing the new 
variable 

4 

r = -£\xi - x 2 \ (11) 
5 

so that 



S 3 (r)=r (12) 

we have 

One can see that the relation (12) is satisfied by this equation. The corresponding equation 
for the probability density P(Au,r) = P(U,r) is: 

The equation for the exponents £ n following from (13) is: 

«l + ^+6,)-f(2 + ^) = (14) 
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Thus, we see that all exponents are determined in terms of £2- Substituting n = 2 
into (14) we derive £ 2 = 0.7257. Using this value we have: 

f n « -1.1815 + (1.3958 + 1.1210 n)5 (15) 

This formula gives the following values for the exponents: £1/4 ~ 0.1131; £1/2 ~ 0.2171; 
£ 3/4 « 0.3140; 6 « 0.4046; £ 2 = 0.7257; £4 = 1.2433; £ 5 « 1.4644; £ 6 « 1.6684. It is 
interesting that, strictly speaking, the scalings of all moments, including those with n < 1, 
are anomalous and cannot be obtained from dimensional considerations. It is possible to 
formally continue the expression (15) into the interval n < 0. We can see that the moments 
S n with n < n c = -fffff ~ -1-2451 do not exist. Since the PDF P(Au, r) ^ at the origin 
Au = 0, it is expected that the moments S n (r) with n < — 1 diverge. Thus, the fact that 
the relation (15) gives the critical moment-order n c —1.24, close to n c — — 1, indicate that 
formula (15) might give close-to- correct values of the scaling exponents £ n for n > n c , which 
are not too far from critical n c = —1. We have: £-1/4 ~ —0.126; £-1/2 ~ —0.269; ^-3/4 ~ 
-0.438; w -0.660 . 

Let us now consider the expression for the dissipation term D which is proportional to 
vq — > 0. The eddy viscosity approximation limits the renormalized perturbation expansion 
in powers of the non-linearity by one-loop contributions only. The technical details of the 
procedure vere described in Refs. [10], [11], [12]. We are interested in the correlation function, 
similar to D: 

D 1 = u <u(q l )....q 2 i u { q i )...u{k - q 1 .... - g n -i) > (16) 

where % — 1; 2; n — 1; there is one contribution involving q 2 n = (k — q\.... — q n -\) 2 . The 

expansion is generated by the iteration procedure involving the Navier-Stokes equations (5) 
or (6) (see Refs. [10], [11]). In the zero -loop approximation we have 

D° w n£S n - 2 (r) 

where the 0(1) dissipation rate £ = u J k 2 E{k)dk. This term is obtained from (16) in the 
limit of the "free momentum" k — > 0. Substituting this expression into the equation of motion 
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(9) we find that, for the solution found above, balance is impossible since — 1 + y + 7^ £n-2- 
Moreover, it is easy to see that the zero- loop contributions to the expression for D, we are 
interested in, cancel in the limit k — > 0. In the one- loop approximation we can have n u.v.- 
divergent contributions of the order v^Snk^ and only one term coming from the /c-dependent 
q n in relation (16) (here A; is a free "momentum"): 

2 ,dS n v[r)S n , 
D i (XV ( r )^~I 0C —^- ( 17 ) 

where the effective (eddy) viscosity u(r) = 0(r 1+ ~). Substituting (17) into (9) and fixing 

the proportionality coefficient to satisfy Kolmogorov's 4/5 law (1) we derive the equation 

(13) used for evaluation of the exponents £„. The u.v.-divergent contributions are assumed 

to cancel for the eddy viscosity approximation to work. This fact was proved to be correct 

for the model case of a simple effective forcing function in Ref. [12]. It was shown that this 

cancellation is the result of overall energy balance that requiers < S ■ u >= £. In general, 

the rigorous demonstration that it is so is a difficult task. All we can say now is that the 

right side of equation (13) is consistent with the eddy viscosity approximation. 

The fact that eddy viscosity works so well for the description of quite complex engineering 

flows was known for many years in the engineering community. For these many years, the 

eddy viscosity concept was treated with suspicion or even contempt by many physisists due 

to the long- standing belief that, being a result of one- loop closure, it cannot be useful for 

the description of the "non-perturbative" intermittency of turbulence. Recently, it has been 

shown in direct numerical experiments on the two-time, two-point correlation function in 

one-dimensional Kolmogorov turbulence for a forced Burgers equation that eddy viscosity 

works well to descibe the small scale dynamics in this complex and strongly intermittent 

system [13]. The work of Shraiman and Siggia [8], as well as other recent contributions 

to the theory of a passive scalar [2]- [4] and [7], also based on the eddy diffusivity, showed 

that anomalous scaling can be derived in the eddy diffusivity approximation, provided the 

zero modes are treated with due respect. (In the problem of a passive scalar advected by a 

white-in-time random velocity the eddy diffusivity is an exact consequence of the equations 

of motion). 
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The scaling exponents, calculated by the present anlysis, are compared with the results 
of physical experiments (Refs. [14]- [16]) in Table I. The measurement of the two-point corre- 
lation functions is very difficult. That is why usually the single-point, two-time correlations 
functions F(t) =< u(x,t)u(x,t + r) > are measured and a simple space-time tranformation 
r = Ut is used for the interpretation of the data in terms of spatial correlation functions. 
This transformation is called Taylor frozen turbulence hypothesis, which is accurate when 
the ratio of the fluctuating u rms and mean U velocities is very large a = U/u rms » 1. 
This criterion is well satisfied in the wall-bounded flows (boundary layers, pipes, ducts etc 
), while in the jets, mixing laers and other open flows a = 0(1). The experimental data 
used in the Table were measured in the atmospheric boundary layer [14] with a m 30, in 
the flow between counterrotating disks [15] with U — and in the turbulent jet [16]. To 
account for the space-time relation in an accurate way, the authors of Ref. [15] introduced 
a Lagrangian-like transformation, based on the idea of the "local Taylor hypothesis", which 
can be accurate even when a = 0. The data of Ref. [16] were collected in the jet flow using 
a novel nonintrusive optical technique which enables one to measure directly multi-point 
spatial correlation functions. Due to the experimental uncertainty, we have avoided com- 
parison with the data on the two-time correlation functions, based on the traditional Taylor 
hypothesis, obtained in the open flows. 

The fact that the scaling exponents evaluated in this work agree so well with experimen- 
tal data is additional evidence that the eddy viscosity approximation is much more powerfull 
than may have expected for a simple one-loop theory. The eddy viscosity approximation, ap- 
plied to the inertial range dynamics, accounts for pressure fluctuations only in the numerical 
value of the factor T in the eddy viscosity definition [10], [12]. This is not good enough for 
the correct representation of the effects ifluenced by large-scale velocity fluctuations, which 
dominate very high-order moments of velocity differences (vortex filaments etc). Moreover, 
the equation (13) takes into account only the first dissipative anomaly given by (1). It is 
clear the anomalies, describing constant or close-to-constant fluxes of other flow character- 
istics, like K = UiUi, can influence the properties of high- order moments. Since the theory 
does not account for these effects, it cannot be used for evaluation of the exponents when 
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n is large enough. To conclude, we would like to mention recent multifractal model by She 
and Leveque [17], leading to a different expression for the scaling exponents £ n which are in 
a good agreement with experimental data for n > 0. The relation between the equation (13) 
and formula (15), derived here, and the She-Leveque theory is not understood. 
I would like to thank K.R.Sreenivasan for communicating his exprerimental findings on the 
scaling exponents of low-order structure functions prior to publication. Helpful discussions 
with R.H. Kraichnan, M. Nelkin and S. Orszag are aknowledged. This work was supported 
in parts by ONR, AFOSR and ARPA. 



n 


£ n ,calc 


a, [14] 


&,[15] 


U16] 


0.1 


0.0466 


0.043 ±0.006 






0.2 


0.0913 


0.083 ±0.010 






0.3 


0.1346 


0.123 ±0.110 






0.5 


0.2171 


0.200 ±0.150 






1 


0.4046 


0.384 ±0.023 


0.40 




1.5 


0.5727 


0.555 ±0.024 






2 


0.7257 1 


0.714 ±0.025 


0.71 


0.70 ±0.01 


4 


1.2433 


1.21 


1.24 


1.28 ±0.03 


5 


1.4644 


1.53 


1.48 


1.50 ±0.05 


6 


1.6684 


1.66 


1.69 


1.75 ±0.10 


8 


2.0378 


2.05 






10 


2.3690 


2.38 







references 

1. A.N.Kolmogorov, Dokl. Akad. Nauk SSSR, 30, 299 (1941) 

2. R.H. Kraichnan, Phys.Rev.Lett., 72, 1016 ( 1994) 

3. K. Gawedzki and A. Kupiainen, Phy. Rev. Lett., 75, 3608 (1995) 

4. M. Chertkov, G. Falkovich, I. Kolokolov and V. Lebedev, Phys.Rev.E, 51, 4924 (1995) 

5. Ya. G. Sinai and V.Yakhot, 1988, unpublished 



9 



6. R. H. Kraichnan, Phys. Fluids, 11, 945 (1968) 

7. V. Yakhot, Phys.Rev.E, 1996 (in press) 

8. B. Shraiman and E. Siggia, C.R.Sci. (Paris) 321, 275 (1995) 

9. R.H. Kraichnan, Phys. Fluids, 9, 1728 (1966) 

10. V. Yakhot and S. A. Orszag, Phys. Rev. Lett., 57, 1722 (1986) 

11. V. Yakhot Orszag, J. Sci. Comp., 1, 3 (1986) 

12. V. Yakhot and L.M. Smith, J.Sci.Comp., 7, 35 (1992) 

13. A. Chekhlov, V. Yakhot, Phys. Rev. E, 51, R2739 (1995) 

14. The low-oredr moments S n < 2 were measured recently by B. Dhruva and K.R. Sreeni- 
vasan, private communication, The moments with n > 2, quoted in the Table, are from G. 
Stolovitzky, K.R. Sreenivasan and Juneja, Phys. Rev. E, 48, R3217 (1993) 

15. J.-F. Pinton and R. Labbe, J.Phys. II France, 4, 1461 (1994) 

16. A.Nullez, U. Frisch, R. Miles, W. Lempert, private communication 

17. Z.-S. She and E. Leveque, Phys.Rev.Lett., 72, 336 (1994) 



10 



